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ABSTRACT 

Aims. We analyse the spectral energy distributions (SEDs) of the most distant galaxies discovered with the Hubble Space telescope and 
from the COSMOS survey and determine their physical properties, such as stellar age and mass, dust attenuation, and star-formation 
rate. 

Methods. We use our SED fitting tool including the effects of nebular emission to analyse three samples of z ~6-8 galaxies with 
observed magnitudes ~ 23 to 29. Our models cover a wide parameter space. 

Results. We find that the physical parameters of most galaxies cover a wide range of acceptable values. Stellar ages, in particular, are 
not strongly constrained, even for objects detected longward of the Balmer break. As already pointed out earlier, the effects of nebular 
lines significantly affect the age determinations of star-forming galaxies at z ~ 6-8. We find no need for stellar populations with 
extreme metallicities or other non-standard assumptions (IMF, escape fraction) to explain the observed properties of faint z-dropout 
galaxies. Albeit with large uncertainties, our fit results show indications of dust attenuation in some of the z ~ 6-8 galaxies, which 
have best-fit values of Ay up to ~ 1. Furthermore, we find a possible trend of increasing dust attenuation with galaxy mass, and a 
relatively large scatter in specific star-formation rates, SFR/Mt,. 

Conclusions. The physical parameters of very high-z galaxies may be more uncertain than indicated by previous studies. Dust atten- 
uation seems also to be present in some z ~ 6-8 galaxies, and may be correlated with galaxy mass, as is also the case for SFR. 
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1. Introduction 

Finding and studying the most distant galaxies formed during 
the epoch of reionisation, more recent than 1 Gyr after the Big 
I Bang, is one of the challenges of contemporary observational 
: astrophysics. Over the past few years considerable progress has 
been made in this field, pushing the observable limits beyond 
redshift 6 with the use of ground-based facilities and satellites. 

A variety of observational programs have tried to locate 
z > 6 galaxies using different observational techniques, mostly 
involving either searches for Lyes' emission through narrow-band 
filters or searches using the Lyman break technique - also called 
the dropout technique. These have been performed either in 

■ blank fields or in fields with galaxy clusters, which act as strong 

■ gravitational lenses, targetting different depths and survey areas. 
The objects found in this way are line emitters or Lyman break 
galaxies (LBGs). 

Although Lya emitters are among the most dista nt galax- 
ies with spectroscopically confirmed redshifts (see live et alJ 
[2006; Otaetal. 2008), few have been fou nd at z ^ 7 (see 
e.g. Cubv et al. 2007; Stark et al. 2007^ IWilhs et alj l200l 
iHibon et al..2009: Sobral et aL.2009.) . Furthermore, because of 
their faintness the photometry available is inadequate in terms 
of depth to allow studies of their stellar populations. 

Surveys using strong gravitational lensing were among 
the firs t to pave th e wa y in the o f stu dy z > 6 galax- 
ies (see 'Kneib et al."2004J, 'Pello et al.H2004!; lEaami et al.''20()5l 
[Richai-d et al. 2006, 2008; Bradlev et al. 2008; Zheng et ay 



a handful of z ~ 7 candidates in blank fields dBouwens et alJ 
12004; Labbe et al..,2006.; .Henrv et al...2008.) . These pilot studies 
also showed that some of the z ^ 7 galaxies could be detected 
at 3.6 and 4.5 yum with Spitzer, thus probing the rest-frame opti- 
cal em ission from these objects (Egami et al...2005. ; Labbe et al] 
120061) . 

Since then, surveys of z-dropout galaxies (targeting z ~ 7 
objects) have been extended to cover larger areas, primarily with 
groun d-based instruments (fMannucci et al. 2007; Capak et aL 
2009; 'Castellano et al.' 2009; Hickev et al. 2009; Ouchi et a" 
2009; Wilkins et al..2009). b ut also with HST (Henrv et al..2007r 
120091: iGonzalez et al.ll2009l) . In most cases, however, only a few 
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12009 ). Ulti-a-deep fields with the Hubble Space Telescope (HST) 
including near-IR observations with NICMOS have uncovered 



near-IR photometric bands are available, providing so far in- 
formation only on source counts and luminosity functions, but 
precluding more detailed studies of the physical properties of 
the so urces. Notable exceptions are the work of ICapak et al.l 
( l2009h . who present three bright ( J ~ 23) z ^ 7 galaxy candi- 
dates from the COSMOS 2 square degree field, and the study of 
Gonzalez et al. (2009) finding 11 fainter (Jiiow~ 26.-27.5)z ~ 7 
galaxies in the two GOODS fields. Both studies benefit from a 
coverage including optical, near-IR, and Spitzer bands. 

Observations taken recently with the newly installed WFC3 
camera on-board HST have just been released, resulting in pub- 
lications from four independent groups identifying faint (Ji25~ 
27-29) z ^ 7 galaxies, based on the combination of the deep- 
est available ACS/HST and WFC3 data (Oeschetal. 2010|; 
'Bouwens et al." 2010a; ' Bunker et all 120091; IMcLure et al.l 120091; 
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Jfan et al. 2009). While these objects are too faint to be de- 
tected at the current limits of t he deepest Spitzer images, a stack 
of 14 z-dropout galaxies from lOesch et al.l (1201 Oi) shows tenta- 
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tive (5.4 and 2.6 c r) detections at 3.6 and 4.5 fim, respectively 
(lLabbeetal.1120101) . 

Given these detected z ~ 1 galaxies (or candidates) with 
available multi-band photometry, it is of interest to determine 
their physical properties such as stellar ages, reddening, stel- 
lar masses, star-formation rates, and related properties such as 
their formation redshift, specific star-formation rate, and oth- 
ers. Several studies have addressed t hese questions using dif- 
feren t modehng tools (see Bouwens et al.l 2010bt ICapak et"an 
l2009t iGonzalez et al.ll2009t iLabbe et al.ll2010l) . However, some 
consider only special types of star-formation histories (con- 
sta nt star-forma tion rate), or zero dust extinction, and except 
for ICapak et"aD (12009.) none of them accounts for the effects 
of nebular emission (hnes and continua) present in star-forming 
galaxies. Neglecting the latter may in particular lead to systemat- 
ically older stellar ages, to low er dust extinction, and differen ces 
in stellar masses, as shown by ISchaerer & de BarrosI (l2009l) for 
z ~ 6 galaxies. Furthermore, the uncertainties in the derived 
physical parameters are not always determined or addressed. 
Last, but not least, no "uniform" study of the entire data sets of 
z ~ 7 galaxies has yet been undertaken using the same method- 
ology and modeling tools. For all these reasons, we present a 
critical analysis of the physical properties of the majority of z ~ 
6-8 galaxies that have been discovered recently. 

Nebular emission can significantly alter the physical param- 
eters of distant star-forming galaxies derived from broad-band 
photometry. The main reason for this is that the emission lines, 
which are invariably present in the H ii regions accompanying 
massive star-formation, strengthen with redshift, because their 
observed equivalent width scales with (1 + z). Since the main 
emission lines are in the optical (rest-frame) domain and few 
are in the UV, their presence can mimick a Balmer break in ab- 
sorption, a signature usuall y interpre ted as an age indicator for 
stellar populations ( Kauffm ann et al.ir2 003: Wiklind et al. 2008). 
This effect of emission lines, and to a lesser extent also nebu- 
lar continuum emission, can lead to degeneraci es in broad-band 
SED fits o f high-z galaxies as e.g., show n by IZackrisson et alJ 
(l2008i) and ISchaerer & de Barrosr (l2009^. The presence of both 
nebular lines and continua and their contribution to broad-band 
photometry is well known in nearby star-forming galaxies, such 
as very metal-poor objects (e.g., I Zw 18, SBS 0335-052, and 
others), blue conipact dwarf galaxies and related objects (cf . 
Izotov et al. 1997'; Pap aderos et"aDl2002t iPustilnik etalll2004t 
Papaderos et al. 2006). The strongest evidence of a significant 
contribution of the nebular continuum in some nearby star- 
forming gala xies is the o bservational finding of a Balmer jump in 
emission (see lGuseva e t al. 2007). For these reasons, it is impor- 
tant to include nebular emission in SED fits of distant starbursts 
and to examine their effect on the derived physical properties. 

In the present paper, we analyse samples of z ~ 6-8 galaxies 
discovered recently. The data are compiled from the literature, 
including the brightest objects from the sample offc apak et al. 
(12009 ). the "intermediate" sample of Gonzalez et al.l (1200 9). and 
the faintest z-dropouts recently found with the WFC3 camera. 
Applying our up-to-date spectral energy distribution (SED) fit- 
ting tool, we search in particular for possible trends in the phys- 
ical parameters of z ~ 7 galaxies over a range of ~ 6 magni- 
tudes, i.e., a range of ~ 250 in flux. First results from our analysis 
are presented here. A more detailed and extensive study of the 
properties of z-dropout galaxies and comparisons with objects at 
lower redshift will be published elsewhere. 

In Sect. 12] we summarise the galaxy sample and the SED 
fitting method. In Sect. [3] we present our results for the three 
subsamples. The overall results of the whole z ~ 7 LBG sample 



and implications are discussed in Sect.|4] where we also compare 
our results to those for LBGs at lower redshift. Our main con- 
clusions are discussed and summarised in Sect. |5] We assume a 
flat ACDM cosmology with Hq = 70 km s"' Mpc ', Qm = 0.3, 
and Qvac = 0.7. All magnitudes are given in the AB system. 



2. Observational data and modelling tools 

2.1. z~l galaxy samples 

To determine the physical properties of z ~ 7 galaxies and their 
uncertainties, we chose the following three samples: 

- Two of the three bright (7 ~ 23) z dropout ga l axies from 
the COSMOS survey, discovered bv lCapak et al] (l2009l) . We 
refer to these as the "bright sample". 

- The 1 1 Z850LP dropout objects identified by iGonzalez et al.l 
(1 20091) from the HST ACS and NICMOS data in the GOODS 
and HUDF fields, plus their mean SED. These objects typi- 
cally have Jiiow~ 26.-27.5, and are referred to as the "inter- 
mediate sample". 

- The "faint sample", including 15 of th e 16 zgjoLp dropout 
candidates found by lOesch et"an (12010.) in the HUDF us- 
ing the newly installed WFC3 camera of HST, and the 
15 additional objects identified as z ~ 6-9 candidates 
by [ McLure et al. (2009). The photometry is taken from 
iMcLure et al.i (l2009l) . They typically span a range from 
Ji2 5~ 27 to 29. We al so include the stacked S ED obtained 
bvlLabbe et all (l2010l) for 14 objects from the lOesch et aP 

sample, which shows tentative (5.4 and 2.6 cr) detec- 
tions in the 3.6 and 4.5 jjm bands of Spitzer. 



The following photometric data/filters was used for the sam- 
ples: (1) /+, z+ from SuprimeCam on SUBARU, J, H, K from 
WIRCA M on the CFH T, and channels 1-4 of IRAC/Spitzer for 
the Capak et al.l (l2009h sample. Since object 2 is detected in the 
i+ band and at 24 nm, and its SED indicates a low redshift 
(z ~ 1.6. ICapak etani2009[) . we exclude it from our analysis. 
(2) B435, Veoe, 1776, zgsoLP filters of ACS/HST, Jnow and Hi6ow 
ofNICMOS/HST,/:^ and channels 1-2 of IRAC/Spitzer for the 
IGonzalez et al.l (l2009l) objects. We adopted the properties of the 
Ks filter of ISAACA^LT for all objects. (3) B435, Veoe, 1776, zgsoL? 
filters of ACS/HST, Y105, J125, Hieo filters of WFC3/HST, and 
channels 1-2 from IRAC/Spitzer for the faint sample. The origi- 
nal photometry from the respective papers was adopted. 

Except for the 3 objects in the bright sample for which one 
spectral line was found for each of them, no spectroscopic red- 
shifts are available for these objects. We therefore treat the red- 
shift as a free parameter for all objects. 

2.2. SED fitting tool 

To analyse the broad-band photometry, we us e a modified ver- 
sion o f the Hyperz photometric redshift code of Bolzonella et al.l 
(I2OOOI) described in ISchaerer & de BarrosI (12009^ The main im- 
provement with respect to both earlier versions and other SED 
fitting codes is the treatment of nebular emission (lines and con- 
tinua), which can have a significant impact on the broad-band 
photometry of high redshift gal axies and hence their derived 
properties (see ISchaerer & de Ba rros 2009). We use a large set 
of spectral templates (primarily the GALAXEV synthesis mod- 
els of Bruzual & Chariot 2003), covering different metallicities 
and a wide range of star formation (SF) histories (bursts, expo- 
nentially decreasing, or constant SF), and we add the effects of 
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nebular emission. Models with a more sophisticated description 
of stellar populations, chemical evolution, dust evolution, and 
different geometries (see e.g.. Schurer et al. 2009.) are not used, 
given the small number of observational constraints. 

We adopt a Salpeter IMF from 0.1 to 100 Mq, and we ac- 
curately consider the returned ISM mass from stars. Nebular 
emission from continuum processes and lines is added to the 
predicted spectra from the G ALAXEV models, as described in 
ISchaerer & de Barrosl (l2009t) . proportional to the Lyman contin- 
uum photon productio n. The relative line intensities of He an d 
metals are taken from Anders & Fri tze-v. AlvenslebenI (l2003h . 
including galaxies grouped into three metallicity intervals cov- 
ering ~ 1/50 Zq to Zq. Hydrogen lines from the Lyman to the 
Brackett series are included with relative intensities given by 
case B. Our treatment therefore covers the main emission lines 
of H, He, C, N, O, and S from the UV (Lya) to 2 jum (rest-frame), 
necessary for fitting the SED of galaxies at z > 4 up to 10 yum 
(IRAC channel 4). 

The free parameters of our SED fits are: the redshift z, the 
metallicity Z (of stars and gas), the SF history described by 
the timescale r (where the SF rate is SFR oc exp"'^'"), the age 
t defined since the onset of star-for mation, the extinctio n Ay 
described here by the Calzetti law dCalzetti et al.l '2000). and 
whether or not nebular emission is included. In some cases, 
we exclude the Lya line from the synthetic spectra, since this 
line may be attenuated by radiation transfer processes inside the 
galaxy or by the intervening intergalactic medium. 

Here we consider z e [0, 12] in steps of 0.1, three metallici- 
ties Z/Zq =1, 1/5, 1/20, T = 5, 7, 10, 30, 50, 70, 100, 300, 500, 
700, and 1000 Myr in addition to bursts and SFR=constant, ages 
up to the Hubble time, and Ay = 0-2 (or 4) mag in steps of 0.2. 
In general, the combination of all parameters leads to ~ 3 x 10'' 
models for each object. 

Non-detections are included in the SED fit with Hyperz by 
setting the flux in the corresponding filter to zero, and the error to 
the Icr upper limit. For all the above combinations we compute 
the and the scaling factor of the template, which provides 
information about the SFR and M*, from the fit to the observed 
SED. Minimisation of over the entire parameter space yields 
the best-fit parameters. 

To illustrate the uncertainties in the resulting fit parameters, 
we examine the distribution of across the entire parameter 
space. To determine confidence intervals from the x^ distribu- 
tion, the degree of freedom must be known to determine the Ax^ 
values corresponding to different confidence levels, or Monte 
Carlo simulations must be carried out. In any case, the pho- 
tometric uncertainties, typically taken from SExtractor, would 
also need to be examined critically, since these may be underes- 
timated, and since errors in the relative photometric calibration 
between different telescopes/instruments, which affect SED fits, 
are usually not taken into account. We chose to plot the ID 
distribution for the parameter of interest, marginalised over all 
other parameters, so that the reader is able to appreciate these 
distributions. Illustrative confidence intervals are determined by 
assuming ls.x^ ~ 1, the value for one degree of freedom. This 
should provide a lower limit to the true uncertainties. More quan- 
titative estimates of the uncertainties will be given in a subse- 
quent publication, which will include the analysis of a larger 
sample of LBGs at different redshift. 



3. Results 

3.1. Photometric redshifts: overview of fhe full sample 

The photometric redshifts Z nhot of the objects from the three sub- 
samples were discussed bv Capa k et al.l (l2009h : iGonzalez et ^ 
(2009), and McLure et al. (2009). Since these authors use dif- 
ferent spectral templates and methods, it is useful to examine 
briefly the redshifts we derive from our SED fits, and their 
dependence om nebular emission. Figure [1] shows the best-fit 
model values for Zphot using either standard templates (i.e., ne- 
glecting nebular emission), or including nebular emission (lines 
and continua), and the latter but neglecting the contribution from 
Lyof. Clearly, the contribution of Lya can lead to higher photo- 
metric redshifts, since it can compensate for the drop of the flux 
shortward of Lya, and hence lead to drop-out at higher Zphot- 
With the prescription used for Lyo- in our models (i.e. maxi- 
mum emission according to Case B recombination) this typically 
leads to Azphot ^ 1. I n some cases, e.g., for 4 objects from the 
Gonzalez et a ll (120091) sample, the shift is larger The reason for 
this large shift is the available filter set, which include zgsoLp and 
J I low for this sample, whereas Yios, a filter that is intermediate 
between zggoLP and J125, is available for the WFC3/HST (faint) 
subsample. By including the nebular continuum and all spec- 
tral lines except Lya (blue symbols) one recovers essentially the 
same photometric redshifts as with standard templates. This is 
expected, since the Lyman break — the main feature determin- 
ing Zphot — can only be strongly affected by Lya. 

We compared our photometric redshifts against objects with 
known spectroscopic redshifts, where possible. For a sample of 
B, y, and /-dropouts from the GOODS fields, we find good 
agreement for the majority of objects using the GOODS-MUSIC 
photometry (Santini et al. 2009). For this sample, spanning ob- 
jects with Zspec ~ 4-6, our results are essentially the same 
with/without nebular emission, and with/without Lya. Since 
Lya emission may be weaker than predicted by the mod- 
els, because of the multiple scattering in the presence of dust 
( Verhamme et al. 2008) and/or because of the intervening IGM, 
we subsequently consider models including all nebular lines ex- 
cept Lya. It must, however, be noted that for objects with strong 
Lya emission the true redshift may be higher than Zphot obtained 
from photometric codes neglecting this line. 

3.2. Bright sample 

Our SED fits for these objects yield results (redshift probabil- 
ity d istributions a r id phy sical parameters) broadly in agreement 
with ICapak et"an (I2009!) . which is unsurprising since these au- 
thors also include nebular lines in their analysis using the Le 
Phare code. Since we include SEDs spanning metallicities from 
Z0 to 1/50 Z0 (in contrast to Z© only), we obtain a wider range 
of acceptable fit parameters. 

For example, we find evidence of significant reddening in 
object 1, with Ay ~ 1- 2.6, in ag reement with the best-fit value 
of Ay -1.2 given bv ICapaket a l. (2009). The corresponding 
X^ distribution is shown in Fig. lA.ll The situation is similar for 
object 3, although for more moderate extinction (Ay ~ 0.6-1.6). 

The corresponding range of ages, SFRs, and stellar masses 
for both objects, obtained with and without nebular emission, 
are illustrated in Figs. lA.2] to lA.4l Approximately (within Ax^ ~ 
1-2.3) object number 1 (3) has thee best-fit model parameters 
t ~ 0-30 (10-200) Myr, log(M*) ~ 10.8-12. (10.6-11.6) Mq, 
and SFR ~ 100-10'' (10-10^) Mq yr^^ We note that a very large 
range in SFR is obtained from the SED fits since both the SF his- 
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Fig. 1. Comparison of best-fit model photometric redshifts for 
the three samples using SEDs without (x-axis) and with neb- 
ular emission (y-axis). Circles, squares, and triangles indicate 
the bright, intermediate, and faint samples, respectively. Black 
symbols show the comparison with all nebular lines, i.e. in- 
cluding also Lya, blue symbols with Lya suppressed. The error 
bars shown here denote the 68% confidence interval derived by 
Hyperz from the redshift probability distribution derived assum- 
ing P(z) oc exp(-x^(z)), where x^(z) is the minimum chi-square 
value over all other parameters. Discussion is given in text. 



tory and the extinction are kept free. The high SFR tail is related 
to solutions with high extinction and very young populations, 
where the UV output per unit SFR remains below the equilib- 
rium value reached after typically ^100 Myr. For comp arison, 
using the standard SFR(UV) calibration of lKennicuti ( 1199 8) and 
assuming z = 8, one obtains SFRobs ~ 240 yr~' without 
extinction correction for object 1, and SFRcorr ~ 6500 yr"' 
adopting Ay - 1 .2 and Calze t ti's att enuation law. 

As noted bv ICapak et al.l (l2009h . object 3 is not detected in 
H, possibly indicating a dip in the flux between J and K (see Fig. 
|2|. If real, this dip could be explained by the 2175 A dust absorp- 
tion feature, as also mentioned by Caoak et al. (2009). Using the 
Galactic extinction law from iSeaton C197 9) provides excellent 
fits with Ay ~ 1.6, as shown in Fig.|2] The possible indication 
for a 2175 A dust absorption feature at such a high redshift is in 
contrast t o evidence so far, sug gesting the absence of this feature 
(see e.g.. lMaioUno et al.ll200^ . 



3.3. Intermediate sample 

The redshift probability distribution of the objects from the 
iGonzalez et a l. (2009) sample is shown in Fig. |3] All objects 
have a best-fit photometric redshift at high z (z ^ 6-6.5). 
However, 6 of the 11 objects (ID 4, 5, 6, 7, 8, 10) also have an 
acceptable fit at low redshift (z ~ 1-2) with a probability compa- 
rable to the high-z solution; the most unreliable objects being ID 
6, 7, and 8. The 3 brightest objects in Jnow (ID 1, 9, 11) favour 
clearly high-z solutions, since they provide the largest "leverage" 
on the Lyman-break between Jnow and the optical data. 

Figure|4]shows an example of an object with both acceptable 
low- or high-redshift solutions of similar quality. The observed 
SED is reproduced well by a low extinction, young, starburst 
at high-z (z = 7.14 here) or by a 4.5 Gyr old stellar popula- 
tions with Ay - 0.4 at z = 1.29. Reducing the probability of a 




4 6 
wavelength [/im] 

Fig. 2. Obs erved (blue p oints) and fitted SEDs (solid lines) of ob- 
ject 3 from lCapak et al.l (2009) using our spectral templates with 
nebular emission. Lya emission has been supressed here. Two 
best-fit models are shown comparing the standard (Calzetti) at- 
tenuation law for Ay =1.2 with the Galactic extinction law of 
[Seaton ( 1979). Using the latter with Ay -1.6 reproduces nicely 
the observed dip in the H-band. The errorbars of the observed 
wavelength indicate the surface of the normalised filter transmis- 
sion curve. Upper limits in flux indicate Icr limits. Red crosses 
show the synthesised flux in the filters. 



low-z interloper would obviously be possible with deeper optical 
photometry, which would place tighter constraints on the Lyman 
break. Deep K-band data (not available for this object) or other 
constraints on the shape of the SED between 2 and 3.6 //m, may 
allow us to distinguish between the two solutions shown here, 
and provide stronger constraints on the possible Balmer break - 
hence the age - of this object. 

Object 1 (UDF-640-1417) from iGonzalez et alj (l2009t) is 
such an example, benefiting from deeper optical imaging and 
Ks data, as shown in Fig.|5] The former leads to a well-defined 
and clearly most probable high-z solution at Zphot ~ 6.7-6.9. 
However, here the observed spectral shape between the rest- 
frame UV and optical range (probed by JHK and 3.6-4.5 //m, 
respectivel y) may imply a degene racy between age and extinc- 
tion. While IGonzalez et alj (l2009h fit this SED with a dust-free 
population of several 100 Myr age (cf. our magenta line), we 
obtain a tighter fit with models including nebular emission for 
a young population plus dust reddening (typically f S 10 Myr, 
Ay ~ 0.8-1.2). The distribution of;^^^ for Ay and other parame- 
ters are shown in the Appendix (Figs. lBn to lB.4l ). 

In total, we find three objects (UDF-640-1417, GNS-zD5, 
HDFN-3654-1216) with best-fit solutions for Ay ^ 0.6-1.2 and 
relatively young age (f $ 10 Myr). Incidentally these are the 
three brightest objects in Jnow, which may suggest a trend of 
extinction with magnitude (cf. below). However, the significance 
of non-zero extinction is not very high, in particular for GNS- 
zD5 and HDFN-3654-1216, whe re the 3.6 A^m flux is a ffected by 
a bright neighbouring source (cf. lGonzalez et alJl2009h . 
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ID 8, Zsun, no nebular 



0.2 




0.05 



redshift z 



Fi g. 3. Redshift probabil ity distribution P(z) exp(-x^iz)) for 
all lGonzalez et al.l ("2009) objects as derived from Hyperz includ- 
ing our standard templates. Coloured lines indicate the mean 
SED (green), and selected objects discussed in the text (GNS- 
zD2, GNS-ZD3, GNS-zD4 in red, blue, magenta). 



Considering the entire parameter space for the whole sam- 
ple (cf. Figs. IB.ll to lBT4l l. we find that age and dust extinction 
of most objects are not well constrained, and could reach from 
few Myr up to the age of the universe at that redshift, and from 
Ay ~ to ^ 1.6 mag for some objects. In particular, the data 
does not allow us to conclude that these galaxies show no sign 
of dust extinction. Furthermore, their age and hence formation 
redshift remains po orly constrained. The same is also true for 
the mean SED from 'Gonzal ez et aP (|2009), which yields results 
compatible with those of the individual objects^, as expected (see 
thick line in the plots). For comparison. iOonzalez et al.i (l2()09h 
find fairly old best-fit stellar-mass-weighted ages ~ 200 - 400 
My£] typically. Both the assumption of SFR=const and A,, - 
lead to the highest age, since both effects minimise the ratio of 
the rest-frame visible/UV light, the main age constraint. 

Allowing for wide range of SF histories, variable extinction, 
and for nebular emission yields, on average, a broader age range 
(between a few Myr and up to the maximum age), room for ex- 
tinction up to Ay ^ 1., stellar masses from 10**'^ to 10'" Mq, 
and SFR ~ 2-100 Mq yr In other words, the properties of the 
galaxies from this sample of intermediate brightn ess (Jiiow~ 26- 
28), ar e clearly more uncertain than indicated bv lGonzalez et alJ 
( |2009|) . who consider only a restricted range of the parameter 
space. 



3.4. Faint sample 

Using the first UDF observations taken with the newly installed 
WFC3 camera onboard HST, four studies have ide ntified ~ 11 
to 20 z ~ 7 galaxy candidates (or zgsoLP drop-outs) dOesch et al.l 
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Fig. 4. Observed (blue poi nts) and fitted SEDs (solid lines) 
of object 8 (GNS-ZD4) from iGonzalez etalJ (l2009l) using stan- 
dard Bruzual & Chariot solar metallicity models. The best-fit 
(X^ = 0.3) is at z = 7.14 with a stellar population of ~ 130 Myr 
(black line), the secondary solution (chi^ - 2.1) is a 4.5 Gyr old 
population at z = 1.29 (magenta). The errorbars of the observed 
wavelength indicate the surface of the normalised filter transmis- 
sion curve. Upper limits in flux indicate Icr limits. Red crosses 
show the synthesised flux in the filters. 



120 lot iBunker et al.ll2009t iMcLure et al.ll2009t lYan et al.ll2009h . 

What can be said about their physical properties? 

3.4.1. Photometric redshifts 



' With their definition, one has f„, = 0.5 x t for SFR=const assumed 
by these authors, where t is our definition of the stellar age. 



iMcLure et al.l (l2009h had previously examined the photomet- 
ric redshifts and uncertainties for their sample, which also 
covers the majority of z ~ 7 galaxies found by the other 
groups dOesch et al.ll2010HBunker et al.ll2009l: lYan e t al. 200^. 
Unsurprisingly, our results usin g a mo dified version of the 
Hyperz code also used by McLure et al.l (12009) and a slightly 
more extended template librar y, confirm th e ir find ings. In par- 
ticular, for the z ~ 1 sample of lOesch et al.l (1201 Oh we find that 
their objects consistently show photometric redshifts with well- 
defined probability distributions peaking between z ~ 6.3 and 
7.6. For fainter z-dropouts and Y-dropouts, the photometric red- 
shift becomes far more uncertain, and a significant fraction of 
the objects could also be low-z galaxies. As already pointed out 
bylMcLure et al. (2009) and Capak et al. (2009), the depth of the 
optical imaging becomes the limiting factor for objects that faint 
in the near-IR. 



3.4.2. UV slope 

One group pointed out that the fainter of these objects had 
very blue UV-continuum slopes, /?, indicative of "non-standard" 
properties of these galaxies (Bouwens et al. 2010b). Their data, 
shown as red squares in Fig.|6] exhibits a trend of decreasing /3 
(as estimated from their (J125-H160) colour) towards fainter mag- 
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Fig. 5. Observed (blue points) an d best-fit SEDs (solid lines) 
of object 1 (UDF-640-1417) from [Gonzalez et alJ (1200^ using 
standard Bruzual & Chariot solar metallicity models. The ma- 
genta line shows the best-fit for templates without nebular emis- 
sion, and assuming that SFR=constant and A,, = following 
[Gonzalez et al. (2009). The age of the population is found to be 
~ 500 Myr. The black line show the best-fit allowing for nebu- 
lar emission and arbitrary SF histories and extinction, yielding a 
much younger age plus some extinction (f = 6 Myr, Ay - 1 .2). 



nitudes. From the very s teep slopes (i.e. , low values of y6 3) 

reached in faint objects, iBouwens et al.l ((2010 b) claim that ex- 
tremely low metallicities and large Lyman continuum escape 
fractions seem to be required to understand these objects, since 
"standard" evolutionary synthesis models predict minimum val- 
ues of y6 ~ -2.5 for young stellar populations. As the data and 
the errorbars from different groups plotted in Fig. [6] show, we 
cannot reach similar conclusions, given the uncertainties in the 
colour measurement used to determine /3. For the bulk of the 
sources, the y6-slope is compatible within Icr with normal val- 
ues of j0 ~ -2.5 or flatter slopes. Furthermore, it is unclear 
whether the (J125-H160) colour exhibits any systematic trend to- 
wards fainter magnitudes. The observations do not exclude dif- 
ferent properties such as extremely low metallicities and large 
Lyman continuum escape fractions for some of the objects at 
z ~ 1 ■ However, the low significance of these deviations do not 
justify making assumptions that differ significantly from those 
commonly adopted for the analysis of lower redshift objects. 

With the WFC3 filters used in this survey, it is possible 
to generate unusually blue (J125-H160) colours in certain cir- 
cumstances from spectral templates including nebular emission. 
Such a case is illustrated in Fig. [7] showing a fit to object 2502 
from the McLure et al. (2009) sample. Shown here is a model 
of a very young stellar population with solar metallicity includ- 
ing lines and nebular continuum emission redshifted to either 
z - 6.96 (dashed, magenta line) or z = 7.97 (solid, black 
line). In the latter case, the strong intrinsic Lya emission (with 
VK(LyQ')'''=" 200 A) boosts the flux in both J 125 and Y105 since 
these filters overlap by » 0.1 fim. This provides a very blue (J125- 
H160) colour, and enough flux in ¥105 to ensure that this object 



Fig. 6. UV slope p between w 1550 and 1940 A computed from 
yS = 4.29(Ji25-Hi6o) - 2.0 for the z ~ 1 galaxy candidates 
from different WFC3/UDF s amples. The photom etry from var- 
ious sources has be en used: lOesch et"al] ( 2010h (red sq uares). 
Bunker et al. ('2009) (yellow triangles). lYan et all (l2009h (green 
circles), McLure et al. ( 2009.) (blue squares, only objects in com- 
mon with Oesch et al.). 



does not appear as a Yios-drop even at z ~ 8. If at z ~ 7, as 
impUed by SED fits excluding the Lya line (magenta line) or 
from simple colou r-criteria designed to select z ~ 7 galaxies (e.g 
lOesch et aL[|2010i) . the (J125-H160) colour is not affected by Lya. 
However, this is not necessarily the case for all z-dropout galax- 
ies since strong Lya emission, if present, may mimic a lower 
redshift (cf. Sect. 13. lb . The likelyhood of this situation remains 
difficult to establish, especially since Lya may be differentially 
affected by dust, and scattered by the IGM. 

3.4.3. Age and reddening 

For the 15 objects in common between the Oesch and McLure 
samples, we find young stellar populations (f ^ 10 Myr) as best- 
fits, and zero extinction, except for objects 688, 835, and 1092, 
with Ay ~ 0.2-0.6. However, as for the objects from the inter- 
mediate sample, the distribution ofx^ is very flat (cf. Figs. lC.l] 
to lC.4b . allowing a wide range of extinctions (Ay ~ 0-1.2), ages 
of f ~ to several 100 Myr, stellar masses from 10^ to few times 
10'* Mo, and SFRs from 0. 1 (or less) to ^ 200 Mq yr^' , for most 
objects in the faint sample. The wide age range, is possible, e.g., 
since the upper limits at 3.6 and 4.5 fim do not provide a strong 
enough constraint on the optical to UV flux of these faint objects. 
Given the rapid evolution with time in the mass/light ratio in- 
volved here (mostly the UV-optical domain), the uncertainty in 
the ages translates into a large spread in stellar masses, as show 
by Fig. IC.3I The wide range of acceptable SFR values is due 
to both age and SF history (parametrised here by the e-folding 
timescale r) being kept free, in contrast e.g., to commonly used 
SFR(UV) calibrations assuming SFR=constant and ages f ^ 100 
Myr 
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Fig. 7. B est-fit SEDs o f objec t 2502 (photometi-y in blue sym- 
bols) of iMcLure et alJ (1200 9^ for solar metallicity models in- 
cluding nebular emission. The black solid line shows the model 
including all lines (Lya in particular) with a best-fit Zphot=7. 97. 
Red crosses indicate the flux in the filters for this model. The 
dashed magenta line shows the best-fit model excluding Lya, 
found at Zphot=6.96. This object has an observed (J125-H160) - 
-0.5 1 colour, corresponding to /3 - -4.2 Note the strong effect 
Lya may have on (J125-H160) and Zphot, although this object is a 
Z85()LP-dropout, thought to be at z ~ 7. 



Fig. 8. Best-fit SEDs to flie SED of the stack of 14 z-dropout 
galaxies from lOesch et"an (2010) measured by Labb e et aP 
62010) . Photometry is shown by blue symbols; including the 1 cr 
limit for Ks. Red crosses indicate the synthetic model fluxes in 
the filters. The magenta line shows the best-fit using solar metal- 
licity templates without nebular emission, yielding a maximum 
age (t ~ 700Myr), Ay - 0, and other parameters also similar to 
iLabbe et al.l (l2010h . The black line shows the best-fit model in- 
cluding nebular emission (and Lya suppressed), with an age of 
/ - 4 Myr, Av - 0.2, and a stellar mass of ~ 5. x 10^ M©, 
more than 1 dex lower than the mass corresponding to the old fit. 



3.4.4. Physical properties from tine stacked SED 

The individual z-dropout candidates of lOesch et al.l ( 1201 0') and 
[McLure et al. (2009) are undet ected in the deep, a vailable 3.6 
and 4.5 fim Spitzer images, but Labbe et al. stacked the 

image s of 14 of the 16 z-dropout galaxies from lOesch et"al] 
(120101) . obtaining 5.4 and 2.6 cr detections in these bands. Fitting 
the SED of this stack, we find that the physical parameters are 
more tightly constrained, as the thick line in Figs. IC.ll to IC.4I 
show. Overall our best-fit values (with or w ithout nebular emis - 
sion) are very similar to those obtained by iLabbe et"al] (1201 Oh . 
bearing the different definitions of stellar ages t and tw in mind. 
Furthermore, the values of the physical parameters derived for 
the majority of the individual objects is compatible with the val- 
ues determined from the stack. 

An uncertainty remains, however, in the age and conse- 
quently also in the stellar mass determination. We first ob- 
tain secondary solutions with Ax^ ~ 1-2 with young ages (~ 
4-5 Myr) and a small extinction {A^, ^ 0.2-0.4) as we con- 
sider all metallicities, both with our without nebular emission. 
Furthermore, if we suppress the Lya line we improve the fit, 
leading to a best-fit at young ages (~ 2-7 Myr). The correspond- 
ing distributions illustrating these results are shown in Figs. 
IC.5l to lC8l and SED fits from these models are shown in Fig. 
[8] This figure clearly illustrates how an apparent Balmer break 
can be explained by an old population (here t ~ 700 Myr) or by 
ne bular emission from a young er population, as already shown 
bv lSchaerer & de BarrosI (l2009l) . A suppression of the Lya line 
is justified since Lya may be attenuated by radiation transfer 



processes inside the galaxy or by the intervening intergalactic 
medium. In any case, neglecting nebular emission is inconsis- 
tent for spectral templates with recent (S 10 Myr) or ongoing 
star formation. Finally, in comparison with the bright and inter- 
mediate samples we may also question why the faintest z ~ 1 
objects should have the oldest stellar populations, whereas the 
ages of brighter objects are compatible with a broad range of 
ages, including young ones. 

The uncertainty in the age also translates into an uncer- 
tainty in stellar mass. Whereas the estimated aver age mass is 
M*_~ (1 - 2) X 10'' M0 for the old population (cf. iLabbe et al.l 
'2OIO), it is more than a factor of 10 lower for young ages (see 
Fig. IC.7b . since nebular emission contributes partly to the rest- 
frame optical domain. The SFR is, however, hardly affected by 
this uncertainty (see Fig. IC.8I 1. since it is more sensitive to the 
rest-frame UV light present in both young and old star-forming 
populations. In consequence, the specific SFR (SFR/M*) could 
be significantly higher than advocated by Labbe et al. (201(ij). 
Before performing spectroscopy for these objects - a currently 
impossible task - to examine if the 3.6 |/m filter is truly affected 
by emission lines as predicted by the model, the present data 
does not allow us to completely rule out one or the other solu- 
tion. 

4. Discussion 

We discuss the effects of varying the model assumptions on the 
physical parameters derived. We then examine the main derived 
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Fig. 9. Comparison of best-fit values for the age t (upper left), ex- 
tinction Ay (upper right), stellar mass (lower left), and star- 
formation rate SFR (lower right) obtained with different model 
assumptions. The x-axis is the value obtained from our "refer- 
ence model" (including nebular emission but no Lya) except 
for the SFR panel; the other values are plotted using different 
colours: models with r > lOMyr (yellow), models with Lya 
(blue), models without nebular emission (red), and models with- 
out nebular emission and SFR= const (green). For the SFR com- 
parison, we use the SFR(UV) value for the x-axis, and black 
symbols for the reference model. Filled circles, square, triangles 
show the bright, intermediate, and faint samples. The dashed line 
is the one-to-one relation. See discussion in text. 



properties and possible correlations among them for the ensem- 
ble of galaxies studied here, and compare any correlations found 
to those of lower redshift galaxies. Finally, we discuss some im- 
pUcations of our results. 

4.1. Effects of varying model assumptions 

Varying the model assumptions in broad-band SED fits affects 
the inferred physical parameters of distant galaxies in wa ys tha t 
have been discussed in several studies, e.g., by Yabe et al. (2009) 
in qui t e some detail (cf. also Papovi ch et al. 2001; Sawicki et al.. 
I2007t iGonzalez et al.l I2OO9I and others). The most relevant 
assumptions are the star-formation histories, metallicity, the 
inclusion of dust extinction, and the adopted extinction law. 
Furthermore, the inclusion of nebular emission and the assump- 
tions ma de to do so also affect the results as shown he re and in 
(cf. also ISchaerer & de Barrosll2009t lYabe et al.ll2009[ and the 
latter for a discussion of the effect of Ha). 

The impact of different model assumptions on the best-fit 
parameters of our sample is shown in Fig. |9] where we plot the 
values from our "reference model" (including nebular emission, 
Lya suppressed, all SF histories, all ages , all extinction values, 
and all metallicities) on the x-axis, and the same from compari- 
son models on the y-axis. The models we consider here for com- 
parison are: 1) models with t > 10 Myr (yellow), 2) models 
with Lya (blue), 3) models without nebular emission (red), and 



Fig. 10. Best-fit value for the age (upper left), extinction (up- 
per right), stellar mass (lower left), and star-formation rate 
(lower right) as a function of the observed Jab magnitude. All 
best-fit values are taken from the models with t > lOMyr 
(model 1). Filled circles, square, triangles show the bright, in- 
termediate, and fai nt samples. Large open circles show the fi t 



values obtained by Capak et a D (120091) . IGonzalez et all (l2009l) . 
and Labbe et"an (l2010l) for the 2 bright objects and from the 
mean/stacked SED, respectively. The dashed line in the SFR- 
Jab panel shows the standard SFR(UV) line assuming z - 1 and 
no extinction for aU objects. See discussion in text. 



4) same as 3), but for SFR=constant only (green). SEDs without 
nebular emission (3) may be taken as an extreme case of models 
with a very large escape fraction (/esc ~ 1) of ionising radia- 
tion from the Lyman continuum, whereas in the other cases we 
implicitly assume no escape (/esc ~ 0), maximising thus nebu- 
lar emission if young massive stars are present. One may theo- 
retically expect an evolution of /esc with redshift, galaxy mass, 
and other properties, although large differences remai n even 
betwe en simulations (see e.g., Gne din et a l. 200^ Wise & CenI 
[20091 ; iRazoumov & Sommer-Larsenli2010l) . Intermediate values 
of /esc may be included in future models. 

Circles, squares, and triangles in Fig. |9] represent objects 
from the bright, intermediate, and faint samples , respectively. 
In Figs. |9] to [T2I we now include all objects of iMcLure et al.l 
(20 0^ in the faint sample, i.e. 15 objects in addition to the 15 of 
Oes ch et aP (l2010l) . to maximise the sample size and because the 
results with the additional sample do not show noticeable differ- 
ences. We do not show the effect of fixing metallicity, since it is 
small compared to the other effects discussed here. We now dis- 
cuss the dependence of the physical parameters on these models 
assumptions one-by-one. 

4.1.1. Stellar mass 



As noted and discussed by lYabe et aP (l2009h and others, the 
stellar mass is the least sensitive parameter to the model 
assumptions, especially when measurements for the rest frame 
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optical domain are available. This is mainly because it is derived 
from the absolute scaling of the overall SED, and that the mass- 
to-light ratio in the optical does not change much with age and 
star-formation history. 

4.1.2. Age 

In c ontrast, age is the most sensitive quantity (cf. lYabe et aP 
12009 ). From Fig.|9]it is clear that models without nebular emis- 
sion (red and green) yield older ages in m ost cases, as already 
discussed above (cf. ISchaerer & de Barrosi,2009 ). By assuming 
constant SFR, we also tend to infer older ages, as is well known, 
but not in all cases (see green symbols). The inclusion (suppres- 
sion) of Lya can also lead to younger (older) ages, as shown by 
the blue symbols. This is the case for some objects because of the 
increase of Zphot (cf. Sect. 13. lb . which in turn demands a younger 
age (steeper UV spectrum) to reproduce the same observed slope 
in the UV. In other words, to some extent there is also a degener- 
acy between age and redshift0 when Lya is taken into account. 
Naively one may have expected the opposite, namely that the 
inclusion of Lya would lead to older age estimates, since for a 
fixed red shift, the contribution of Lya produces a bluer UV spec- 
trum (cf. 'Finkelst ein et ani2009h . This shows the importance of 
consistently fitting the physical parameters and redshift using the 
same spectral templates. 

4.1.3. Extinction 

Although a spread is obtained in Ay for the different model as- 
sumptions, the best-fit model values of Ay correlate reasonably 
well around our "reference value" for most objects. The greatest 
differences arise when the SF history is varied - here impos- 
ing T > 10 Myr leads to lower Ay and higher ages (because of 
the age-extinction degeneracy), as shown by the yellow sym- 
bols - and, for the same reason, when nebular lines are omitted, 
as shown by some outliers with red symbols. The largest differ- 
ences in stellar mass are found for some objects where nebular 
emission leads to a strong age reduction and hence a lower M/L 
ratio. 

4.1.4. SFR 

The star-formation rate SFR deserves special comment. The 
value of the current SFR (SFR(f)) obtained from SED fits de- 
pends strongly on both SFH and age. Formally, SFR=0 for in- 
stantaneous burst models, which are also considered here and in 
other publications. Furthermore since the UV luminosity emit- 
ted per unit SFR varies by ~ 1 order of magnitude within ~ 
100 Myr even for constant star formation, the cuiTent SFR re- 
sulting from SED fits can diffefl from the usual SFR(UV) cal- 
ibrations by 1 dex or more depending on both SF history and 
age. For this reason, we compare in Fig. |9]the best-fit SFR val- 
ues to_SFR(UV)j[he value derived from the J magnitude using 
the iKennicutl (Il998h caUbration, and assuming no extinction. 
The following differences can be seen: SFR > SFR(UV) is ob- 
tained for many objects, where Ay > and f « 100 Myr. Lower 
SFR(r) <SFR(UV) values are obtained for some objects with 
relative short timescales t, i.e., rapidly declining SF histories. 



^ In addition, there is the well-known age-extinction degeneracy. 

^ In general, one has SFR(0 >SFR(UV), since stars over a broad age 
range contribute to the UV output, and since Lyy/SFR is lower at young 
ages. 



The comparisons shown here all include a marginalisation 
over the three metallicities Z considered in our models. As al- 
ready mentioned, fixing Z leads to small differences. Given our 
limited knowledge of high-z galaxies, other variations in the in- 
put physics, e.g., for the star-formation history or the reddening 
law, could be considered. Quantifying this is beyond the scope of 
this paper. The effect of different reddening laws on the physical 
parameters, howev er, can be un derstood quite simply, and is e.g., 
illustrated in Yabe et al. I (l2009l) . Rising SF histories , which are 
not co nsidered here, have, e.g., been advocated bv iFinlator et al.l 
(l200l . We note, that these SF histories would in general cor- 
respond to a higher extinction th an decreasing on es, to repro- 
duce the observed SEDs of LBGs iFinlator et all (cf. l2007D . This 
would strengthen our arguments in favour of dust at z ~ 7. 

4.2. Correlations between physical properties inz~l 
galaxies 

Although the uncertainties in the physical parameters are rela- 
tively large for most objects, it is helpful to search for corre- 
lations between them, and with observed quantities. Figure [TO] 
shows the best-fit values for age. Ay, stellar mass, and SFR of 
all objects from the 3 samples, plotted as a function of the J- 
band magnitude, which traces approximately the rest-frame UV 
at 1500 A. The values plotted here are taken from our model 1, 
i.e., the reference model plus the constraint that t > 10 Myr, 
in particular to assure that SFR is pro perly defined. We also 
com pare the bes t -fit va lues obtained by Gonzalez et al.l (l2009l) 
and iLabbe et al.l (l2010l) for the mean/stacked SED for the in- 
termediate and fain t samples, and the available fit results from 
ICapak et"an (l2009h for the bright objects. The origin of the dif- 
ferences with these authors have already been discussed above. 

As is clear for the different samples, stellar ages and extinc- 
tion cover a wide range of values with no clear or strong trend. It 
is possible that a tendency of increasing Ay for brighter objects 
exists, as would be expected from other studies at lower redshift 
(cf. below). This trend becomes clearer when Ay is plotted as a 
function of stellar mass, as shown in Fig.[TT] We also find a ten- 
tative trend of the extinction A v with galaxy mass M*, as shown 
in Fig. [n] 

Both the stellar mass and SFR show clear trends with the J- 
band magnitude, albeit with a significant scatter. We note that 
the scatter in the mass-magnitude relation for the z ~ 1 ob- 
jects does not significantly decrease when 3.6 or 4.5 fim pho- 
tometry is used. Deviations of the best-fit SFR from the simple 
"standard" calibration, also indicated in Fig. [10] for A i/ - 0, are 
caused by non-zero extinction, age effects, and exponentially de- 
creasing SF histories, or combinations thereof, as already dis- 
cussed above. Figure [12] shows the corresponding correlation 
between stellar mass and SFR, suggesting a fairly well-defined 
mass-SFR relation at z ~ 7. Our best-fit values yield on average 
higher specific st ar-formation ra t es (SS FR -SFR/M,^) than the 
values derived bv lGonzalez et al.l (l2009l) and lLabbe et al.1 (1201 Oh 
shown by large open circles, and tend to indicate a relatively 
large spread in SSFR. 

4.3. Comparison with lower redshift galaxies 

Although an exhaustive comparison with studies of the physical 
properties of LBGs at lower redshift is clearly beyond the scope 
of the present paper, it is worthwhile comparing briefly our re- 
sults with those obtained at z ~ 5, and at lower redshift. To do so 
we refer to the detailed SED fitting analysis and discussions of 
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Fig. 11. Black circles, green squares, blue triangles show best-fit 
values of Av versus derived from our reference model for 
ojects from the bright, intermediate, and faint samples respec- 
tively. Red filled symbols show the same but from the reference 
model imposing the restriction of r > 10 Myr (model 1). The 
dotted line, given by Ay - log(Mv,/10'*)"'^, is shown to guide the 
eye. Small crosses show the fit results from lYabe et alj d2009l) at 
z~5. 



lYabe et al.l (I2OO9I and references therein), who have compared 
the physical parameters of z = 2, 3, 5, and 6 LBGs from different 
samples. 

The analysis of both Yabe et al. and I Verma et al.l (l2007h of 
z ~ 5 galaxies find clear signs of significant extinction. The for- 
mer study finds E{B - V) values ranging from to ~ 0.5, with a 
medianof£'(B-y)- 0.22, corresponding to Ay - 0.89 i.e. a fac- 
tor 9 attenuation of the UV flux! Our typical Ay values are lower 
and the uncertainties are large; but in many cases the best-fit 
model extinction is Ay > for the z ~ 7 LBGs. For comparison, 
^ouwens et al. (2009, 2010b) advocate an attenuation of the UV 
flux by a factor 1.35-1.6 (corresponding to Ay ~ 0.12-0.19 for 
the Calzetti law) at z ~ 1, and an at tenuation by 2.7 ( Ay ~ 0.4) 
atz~5. The extinction obtained bv lYabe et alj ( l2009h and from 
our analysis is thus typically higher by a factor 3 than in the work 
of Bouwens and collaborators, who estimate the dust attenuation 
from the UV slope. 

We find a trend of increasing extinction with galaxy mass 
(Fig. [TTI) . At z ~ 0-2, it is known that dust extinction in- 
creases with bolometric luminosity (i.e., also with SFR), which 



in turn i ncreases with st ellar mass (e.g..lBuat et al.ll2005 



2008 



20061 



Burgarel la et al.1 120071: iDaddi et al.1 120071: iReddvetal 
2OO81) . Luminosity-dependent dust corrections have been pro- 
posed, e.g., by Reddv & Steidel (2009) and Bouwens et al. 
(I2009h . based on observed variations of the UV slope with Muy. 
From these results, it is not surprising to find a similar rela- 
tion at z ~ 7, here expressed as extinction versus stellar mass. 
The origin of a mass-extinction relation is most likely related 
to that o f the mass-metallicity relation (cf. Tremonti et al. 2004; 
lErb et al. 2006; Finlator et aLl l2007l; iMaioUno et £11120081) . The 
data of lYabe et aL (.2009.) at z ~ 5, also plotted in Fig.[TTl may 



Fig. 12. Stellar mass-SFR relation for the z ~ 1 galaxies. 
Open symbols (circles, squares, triangles) show the "standard" 
SFR(UV) value (not corrected for extinction) versus mass de- 
rived from our reference model for objects from the bright, inter- 
mediate, and faint samples respectively. Red filled symbols show 
the best-fit model SFR and values when assuming t > 10 
My r (model 1). Filled he xago ns denote the mean values derived 
bv iGonzalez et al.l(l2009h andlLabbe et al.1 (120101) . Small crosses 
show the fit results from Yabe et al ■I(l2009h at z ~ 5. T he dashed 
line shows the relation found by Daddi et aTl (|2007|) for z ~ 2 
SF galaxies. The dotted (dash-dotted) lines show the locus for 
SFR=const from z - °° (10) to 7 corresponding to SSFR=1.3 
(3.6) Gyr-i. 



show a less clear trend with mass and some offset. A systematic 
differential analysis of LBG samples at different redshifts, will 
be necessary to shed more light on these issues. 

The best-fit model values of the other physical parame- 
ters (age, mass, a nd SFR) span a similar range as found by 
Yab e et all (l2009h for z ~ 5 L BGs. In coii t rast to the re- 
sults of iGonzalez et alJ (l2009l) and lLabbe etal] (1201 Oh. our rel- 
atively y oung ag es resemble those found bv I Verma et al.l (l2007h 
and Ya be et al.l (|2009) for z ~ 5, whi c h are younger than 
z ~ 2-3 LBGs (c f. Shat^l ev et al.1 12001I llool ISawickietal 



Yabe et aLl2009). However, other stud ies fe.g. JStark et al 



'20071 

2007at lEvles et al.ll2007t fStark et al.ll2009l) . find relatively old 
ages at z ~ 6. and I Yabe et al.l (l2009r confirm some of them with 
their method. Whether LBGs show a clear trend of decreasing 
age with increasing redshift, as one may naively expect, thus re- 
mains to be clarified. 

For the IRAC-detected objects at z ~ 7 (the bright and inter- 
mediate samples), the relation between stellar mass and absolute 
optical magnitude (derived from the 4.5 //m flux) is very similar 
to that of the z ~ 5 objects. The lower masses of the faint z ~ 7 
subsample - undetected by IRAC - provide a natural extension 
towards fainter objects. A similar behaviour is also fo und for the 
SFR, when comparing our z ~ 7 results to those of lYabe et al.l 
(2009). 

The M*-SFR relation found at z ~ 7 resembles that of z ~ 5 
LBGs (see Fig.fTSTi. We note that most objects lie above the M*- 
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SFR relation found at z 2 (cf. lDaddi et alJl2007tlSawicki et alJ 
I2007h . However, both the behaviour and the scatter found for 
the z ~ 7 objects resembles that of z » 5 LBGs, also shown in 
this plot. Given the large uncertainties for the z ~ 7 objects, and 
the different methods used in these analysis, it is possible that 
the same M*-SFR relation (with a similar scatter) is in place 
from z ~ 7 to 2 and that the specific SFR declines t o lower 
redshift (cf. iDaddi et al. 20 07; Elb az et al.. .2007 ; Noeske et alj 
l2007tlGonzalez et alJl2009h . 

A detailed study of the physical parameters of LBGs over 
a wide redshift domain with our modeling tools, including 
more in-depth comparisons of the results from different meth- 
ods/groups, will be presented elsewhere. 

4.4. Implications 

An important result of this study is the possible discovery of dust 
extinction in LBGs at z ~ 7 from detailed SED modeling, which 
contrasts with claims of basically dust-free objects from stud- 
ies of their UV slope or from too restrictive SEP modeling (cf . 
iBouwens et alll2010bt iGonzalez et aLll2009l: iLabbe et al.ll2010l) . 
If this were correct, the UV attenuation - up to Ay ~ 1.2 (factor 
20) for the brightest objects - would imply a higher SFR density 
at high redshift than previously inferred. However, to quantify 
the average dust correction remains difficult, especially given the 
large uncertainties for the individual objects and the possible de- 
pendence of Ay on the UV luminosity and/or stellar mass. Even 
if our results imply that a non-negligible dust correction is re- 
quired at z ~ 7, it is possible that the SFR density decreases with 
redshift from z ^ 3-5 on, at l east if the median attenuation of 
Ay ~ 0.9 found bv lYabe et al.l (l2009l) at z ~ 5 is representative. 

Our analysis, covering a wider domain in parameters space 
than previous studies and also allowing for the effects of neb- 
ular emission, shows that a wide range of stellar ages is ac- 
ceptable for most z ~ 7 LBGs. The Balmer break observed 
for some objects does not always imply old ages, or corre- 
spondingly high formation r edshifts as already pointed out by 
ISchaerer & de BarrosI (l2009h . It can often equally well be ex- 
plained by nebular emission, younger populations at non-solar 
metallicity, extinction, or a combination thereof. 

Younger ages and variations in the star-formation history 
can lead to solutions with lower masses and higher SFR for 

many of the z ~ 7 LBGs than estimated by other groups 

(IGon zalez et al. 2009; Labbe et al. 2010). The specific star- 
formation rate (SSFR=SFR/M*) of these high-z galaxies may 
also be higher than thought. 

In short, our SED analysis including in particular the effects 
of nebular emission and considering a wide parameter space (SF 
histories, extinction etc., shows that the physical parameters of 
z ~ 7 galaxies may differ significantly from those advocated by 
other groups, with possibly many implications for our picture of 
galaxy formation and evolution in the early universe. However, 
we emphasize that the physical properties of these objects are 
affected by large uncertainties (cf. Sect. [3] and Appendix), given 
their faintness and the available data. Furthermore, comparisons 
with results for lower redshifts may be problematic because of 
the use of different methods. 

Obviously, additional and higher quality data, and a detailed 
differential study of the physical parameters of LBGs over a wide 
redshift domain using state-of-the-art tools will be very helpful 
in providing a clearer and more accurate picture of the properties 
and evolution of high-z galaxies in the distant universe, and their 
link to lower redshift galaxy populations. 



5. Conclusions 

We have presented a homogeneous and detailed analysis of 
broad-band photometry for three samples of z ~6-8 galaxies dis- 
covered by the COSMOS survey and with HST ( see Capak et al. 
2009; Gonza lez et alj 120091: iOesch eTdl [20T0I; iMcLure et al. 
,2009; Bunker et alJl2009h . Their 7-band magnitude span a range 
from ~ 23 to 29, the bulk of them having Jab ~ 26-29. 
The broad-band SEDs have been fitted using our modified 
version of the Hyper z pho tometric redshift code described in 
ISchaerer & de BarrosI (l2009i) . which accounts for the effects of 
nebular emission. 

In contrast to earher studies that assumed, e.g., constant star 
formation and/or no extinction, we have explored a wide range 
of the parameter space without using priors. The free parame- 
ters (and range) of our SED fits are: the redshift z ([0, 12]), the 
metalUcity Z (Z/Zq =1, 1/5, 1/20), the SF history described by 
the e-folding timescale t ([0, 1000] Myr, 00), the age f defined 
since the onset of star-formation (< ?//), the extinction A v ([0, 2] 
in general) described by the Calzetti law (Cal zetti et al.ii200Cih . 
and whether or not nebular emission is included^ 

Our main results can be summarised as follows: 

- Overall, we find that the physical parameters of most galax- 
ies studied here cover a wide range of acceptable values, e.g., 
within Ax^ ^ 1 from the best-fit modeling. This finding is in- 
dependent of whether nebular emission is included or not. 

- Stellar ages, in particular, are not tightly constrained, even 
for objects detected with Spitzer, i.e., with photometry both 
blue- and redward of the Balmer break. When nebular emis- 
sion is taken into account, we find that the majority of the 
objects (and the stacked SEDs as well) are most accurately 
reproduced by ages t < 100, which ar e younger than derived 
in other studies o f the same objects ( IGonzalez et al.l 120091; 
ILabbe etani2010h . The younger ages are due to the contri- 
bution of nebular lines to the broad-band rest-frame optical 
filters, which mimic to some extent a Balmer break, as al- 
ready shown bv lSchaerer & de BarrosI (l2009l) . 

- Examination of the UV slope and SEDs of faint z-dropout 
galaxies found with WFC3/HST shows no need for "un- 
usual" stellar populations, extreme metalli cities, or other 
physica l processes, advocated previously by iBouwens et al.l 
(l2010bl). when the uncerta inties are taken into account (see 
also' Finkelstein et all200 9). 

- Albeit with large uncertainties, our fit results show indica- 
tions of dust attenuation in some of the z ~ 6-8 galaxies, 
with best-fit model values of A 1/ up to ~ 1, even among rela- 
tively faint objects (Jab ~ Ib-Tl; cf. Fig.lTOli. 

- We find a possible trend of increasing dust attenuation with 
the stellar mass of the galaxy (Fig.lTTIi and a relatively large 
scatter in specific star-formation rates, SFR/M*. 

- Our results, including the evidence of dust in z ~ 6-8 galax- 
ies, are consistent with the results and tre nds from other SED 
studies o f LBG samples at z ~ 5 (see e.g. lVerma et alj|2007l ; 
lYabe et aL.2009) . 

We will present elsewhere a systematic study of the evolu- 
tion of the physical parameters of LBGs at different redshifts, 
adopting a homogeneous method and including a detailed error 
analysis. 
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since this line may be attenuated by radiation transfer processes inside 
the galaxy or by the intervening intergalactic medium. 



12 



Daniel Schaerer and Stephane de Barros: On the physical properties of z ~ 6-8 galaxies 



Roser Pello for regular discussions and support with Hyperz. This work was sup- 
ported by the Swiss National Science Foundation. 



Capak et al. 2009: ID 1,3 




12 3 



Av 

Fig. A.l. Distribution of the minimum value as a function of 
the attenuation Ay for object 1 (blac k lines) and 3 (blue lines) 
from the bright sample dCapak et al.ll2b09i) . Solid lines show;^^^ 
obtained with standard spectral templates (neglecting nebular 
emission), dashed lines for templates including nebular emis- 
sion, but with Lya suppressed. 
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Fig. A.2. Same as Fig. lA. li as a function of stellar age. 
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Fig. A.3. Same as Fig. lA. li as a function of stellar mass. 



Fig.B.l. Distribution of the minimum value as a function of 
th e attenuation Ay for all objects from the "intermediate sample" 
of iGonzalez et alJ (l2009l) computed with the spectral templates 
including nebular emission, but with Lya suppressed.. Different 
colours represent different objects. The thick purple line show 
for the mean SED. 
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Fig.A.4. Same as Fig. lA.ll as a function of the star-formation 
rate. 



Appendix A: Bright sample jCapak et al .1120091) 

Appendix B: Intermediate sample ^Gonzalez et al. 
12009.) 

Appendix C: Faint sample 

C. 1. \Oesch et all t2010^) objects 
C.2. Stacked SED tLabbe et al\\20ldi) 
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